Experimental evidence of high-resolution ghost imaging and ghost diffraction with 

classical thermal light 
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High-resolution ghost image and ghost diffraction experiments are performed by using a single 
source of thermal-like speckle light divided by a beam splitter. Passing from the image to the 
diffraction result solely relies on changing the optical setup in the reference arm, while leaving 
untouched the object arm. The product of spatial resolutions of the ghost image and ghost diffraction 
experiments is shown to overcome a limit which was formerly thought to be achievable only with 
entangled photons. 

PACS numbers: 42.50.Dv, 42.50-p, 42.50.Ar 



The ghost imaging protocol, provides a flexible way of 
performing coherent imaging by using spatially incoher- 
ent light. It relies on the use of two spatially correlated 
beams, one of them illuminating an object, while the 
other holds a known reference optical setup. Information 
on the object spatial distribution is obtained by correlat- 
ing the spatial distributions of the two beams. Tradition- 
ally parametric down-conversion (PDC) is the source of 
the correlated beams. In the low-gain regime single pairs 
of entangled signal-idler photon can be resolved, and the 
information is extracted from coincidence measurements 
while in the high-gain regime several photon 
pairs form entangled beams and the information is con- 
tained in the signal-idler intensity correlation 0, S • 
Landmark experiments in the low-gain regime showed 
that using entangled photons both the object image could 
be retrived (the "ghost image" experiment 1]) as well as 
the object diffraction pattern (the "ghost diffraction" ex- 
periments 05 0)- 

Recently, a very lively debate arised, aiming to iden- 
tify which aspects (if any) of ghost imaging truly re- 
quires entaglement and which could be reproduced by 
using classically correlated sources. The first theoretical 
interpretation of the experiments suggested that entan- 
glement of photon pairs was essential to retrieve infor- 
mation from the correlations 0, 01 • This interpretation 
was challenged both by theoretical arguments 6] and 
experiments 0, 0, which showed that virtually any 
single result of ghost imaging could be reproduced by 
using classical sources with the proper kind of spatial 
correlation. However, theoretical investigations showed 
that only quantum entanglement provides perfect corre- 
lations in both photon position (near-field) and momen- 
tum (far-field) and suggested that this is crucial 
if both the image and the diffraction pattern are to be re- 
trieved from the same source and leaving the object arm 
untouched lof . Along the same line, recent experimen- 
tal works 0, 0, ^3 , pointed out a momentum-position 
realization of the EPR paradox using entangled photon 
pairs produced by PDC. Based on this result, the Au- 
thors of 0, Q argued that in ghost imaging schemes 



entangled photons allows to achieve a better spatial res- 
olution than any classically correlated beanms, and in 
particular set a lower limit to the product of the resolu- 
tions of the images formed in the near and in the far-field 
of a given, single, classical source. 

This claim is quite in contrast with the theoretical anal- 
yses of , where some of us proposed a particular source 
of classically correlated beams, that can emulate the be- 
haviour of entangled beams in all the relevant aspects 
of ghost imaging, including the resolution capabilities. 
The proposed scheme relied on the use of two beams ob- 
tained by dividing on a beam splitter an intense beam 
with a thermal-like statistics. The two outcoming beams 
were shown to have strong, although completely classi- 
cal, spatial correlations simultaneously in the near and 
in the far-field planes. Thanks to this, the scheme is able 
to reproduce all the results of ghost imaging with entan- 
gled beams, and numerical simulations confirmed this. 
In this paper we provide the first - to our knowledge - 
experimental evidence of the theoretical results obtained 
there. We show that both the ghost image and the ghost 
diffraction results emerge from the correlation measure- 
ments, that we may pass from one to the other by solely 
operating on the setup in the reference arm, and that the 
product of the resolutions clearly overcomes the limit set 
in 0,^^.This definitively demonstrates that entangle- 
ment is not necessary for ghost imaging. 

The experimental setup is sketched in FigQ The 
source of thermal light is provided by a slowly rotat- 
ing ground glass placed in front of a scattering cell con- 
taining a higly turbid solution of 3 /im latex spheres. 
When this is shined with a large collimated laser beam 
(A = 0.6328 /im, diameter Dq ^ 10 mm), the stochastic 
interference of the waves emerging from the source pro- 
duces a time-dependent speckle pattern characterized by 
a correlation time Tcoh on the order of Is. The latter is 
associated to particle motion and can be tuned by vary- 
ing the turbidity of the solution. A small portion of the 
speckle pattern is selected by a I) = 3 mm diaphragm at 
a distance zq = 400 mm from the thermal source. The 
diaphragm realizes an angular selection of the pattern 
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FIG. 1: Scheme of the setup of the experiment (see text for 
details). 



allowing the formation of an almost collimated speckle 
beam characterized by chaotic statistics 15j and by speck- 
les of size Ax ^ Xzo /Dq « 25 fim . The speckle 
beam is separated by the beam splitter (BS) into "twin" 
speckle beams, that exhibit a high (although classical) 
level of spatial correlation. The two beams emerging from 
the BS have slightly non-collinear propagation directions, 
and illuminate two different non overlapping portions of 
the charged-coupled-device (CCD) camera. The data are 
acquired with an exposure time (1 ms) much shorter than 
Tcoh^ allowing the recording of high contrast speckle pat- 
terns. The frames are grabbed at a rate of 1 Hz, so that 
each data acquisition corresponds to uncorrelated speckle 
patterns. 

The optical setup of the object arm 1 is fixed. An ob- 
ject, consisting of a thin needle of 160 /im diameter in- 
side a rectangular aperture 690 /im wide, is placed in 
this arm at a distance di from the BS. A single lens of 
focal length F = 80 mm is placed after the object, at a 
distance pi from the object and qi = F from the CCD. 
In this way the CCD images the far-field plane with re- 
spect to the object. However, the light being incoherent, 
the diffraction pattern of the object is not visible (nor 
its image, in this setup), in the intensity distribution of 
the object beam on the CCD, as shown by FigEh- We 
consider two different setups for the reference arm 2. In 
the first one, an additional lens of focal length F^ is in- 
serted in the reference arm immediately before the lens 
F. The total focal length F2 of the two- lens system is 
smaller than its distance from the CCD q2 = F, being 
jr ~ j; -\- j^. This allowed us to locate the position of 
the plane conjugate to the detection plane, by temporar- 
ily inserting the object in the reference arm and illumi- 
nating it with the laser light. After having determined 
the position of the object that produced a well focussed 
image on the CCD (image shown in FiglSb), the object 
was translated in the object arm. The distances in the 
reference arm approximately obey a thin lens equation of 
the form l/{p2 - di) ^ l/q2 ^ I/F2 The data of the 
intensity distribution of the reference arm are acquired 
and each pixel is correlated with the total photon counts 
of the object arm (i.e. with the sum of the photocounts 
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FIG. 2: Recontruction of the object image via correlation 
measurements (FigQ] with the lens F^ inserted), (a) The 
intensity distribution of the reference arm is cross-correlated 
with the total photocounts in the object arm (statistics over 
5000 CCD frames); (b) image observed by shining the laser 
light on the object, (d) The averages of 500 horizontal sections 
of the images shown in (a) (circles), and (b) (full line) 



over a proper set of the pixels in the object arm), which 
corresponds to having a "bucket" detector in arm 1. Av- 
erages performed over few thousands of data acquisitions 
are enough to show a well-resolved reproduction of the 
image of the needle. This is shown in Fig|2^ and can be 
compared with the image obtained by illuminating the 
object with the laser light (Figl^b). Fig|2t plots the hor- 
izontal sections of the needle image as obtained via cor- 
relation (circles) and with the laser illumination (light 
full line), after averaging over 500 pixels in the vertical 
direction(exploiting the vertical symmetry of the image). 
The spatial resolution shown by correlated imaging with 
incoherent light is comparable with that obtained via co- 
herent illumination. 

In the second setup, without changing anything in the ob- 
ject arm, the lens F^ is removed from the scheme of FigQ 
so that the reference beam passes through a single lens of 
focal F, placed at a distance q2 = F from the CCD. The 
spatial cross-correlation function of the reference and test 
arm intensity distributions is calculated by making aver- 
ages over few hundreds of independent data acquisitions, 
and shows a sharp reproduction of the diffraction pat- 
tern of the object (FiglHt)). This is comparable with the 
diffraction pattern obtained by illuminating the object 
with the laser light (FiglSh) (hence, by simply removing 
the scattering media and observing the light distribution 
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FIG. 3: Recontruction of the diffraction pattern via corre- 
lation measurements (FigQ with the lens F' removed), (a) 
Single-shot intensity distribution in the object arm, showing 
no diffraction pattern, (b) Cross-correlation of the object- 
reference intensity distributions as a function of X2—X1 (statis- 
tics over 500 frames); (c) diffraction pattern observed by shin- 
ing the laser light on the object, (d) Horizontal sections of the 
diffraction patterns shown in (b)(dotted)plotted as a function 
of X = X2 — xi, and (c) (full line) 



of the object arm on the CCD). Horizontal sections of 
these diffraction patterns are plotted in Fig. and 
shows basically that the spatial correlation measurement 
of the speckle beams offers a perfect reproduction of the 
diffraction pattern of this object. Notice that at differ- 
ence from the ghost image experiment a spatial average 
is also involved, which improves the covergence rate as 
described in detail in 0. 

Thus, in this way a high-resolution reconstruction of both 
the image and the diffraction pattern of the object have 
been performed by only operating on the optical setup of 
the reference arm and by using a single classical source. 

The basic theory behind the setup shown in Fig. ^ has 
been explained in detail in Ref. 7] . The starting point is 
the input-output relations of a beam splitter 

bi{x) = ta{x) + rv{x) , b2{x) = ra{x) + tv{x) (1) 

where bi and 62 are the object and reference beams 
emerging from the BS, t and r are the transmission and 
reflection coefficients of the BS, a is the speckle field and 
is a vacuum field uncorrelated from a. The state of a 
is a thermal mixture, characterized by a Gaussian field 
statistics, in which any correlation function of arbitrary 
order is expressed via the second order correlation func- 



r(f,fO = {a\x)a{x')) 



(2) 



The object information is extracted by measuring 
the spatial correlation function of the intensities 
(^1(^1)^2(^2)), where liixi) are operators associated to 
the number of photocounts over the CCD pixel located 
at Xi in the i-th beam, and from this calculating the cor- 
relation function of intensity fluctuations 

G(fi,f2) = - {Il{xl)){l2{X2)) . (3) 

The main result obtained in 7] was 

G{xi,X2) = \tr\'^ J dx[ J dx2hl{xi,x[)h2{x2,X2)T{x,x^) 

where hi and /12 are the impulse response function de- 
scribing the optical setups in the two arms. Let us con- 
sider the scheme shown in FigQ The beams bi and 62 
follow identical paths from the BS to the object plane, 
which we shall refer to as the near-field plane. Due to 
the linearity of the BS transformation, the two beams at 
the near-field plane can be expressed via the same trans- 
formation as (^), with the a replaced by the speckle field 
which has propagated from the source to the near-field. 
In other words, our setup is equivalent to a setup where 
the physical BS is replaced by a point-like beam splitter 
performing the tranformation Q immediately before the 
object. 

The setup of arm 1 is fixed, and by denoting with 
T{x) the object transmission function, hi{xi^x[) = 
{iXF)~^e~^^^'^^T{x[). In the ghost image setup (FigQ] 
with the lens F' inserted), the detection plane of the ref- 
erence arm is the conjugate plane with respect to the 
near-field plane. Hence, apart from inessential phase fac- 
tors h2{x2,X2) = m6{mx2 + ^2), where m ^ 1.2 is the 
magnification factor of the two-lens system. Inserting 
this in Eq. 



G(xi, X2) (X 



J dx^Ti^i^Xi., 



-mf2)T*(f;)e'Xf^i-^i 



\T{-mx2) 



(5) 

2 

(6) 



where r^(x, x^) is the second order field correlation func- 
tion of the speckle beam in the near-field. The second 
line of Eq. Q was derived under the assumption that the 
smallest scale over which the object changes is larger than 
the length over which Tn{x — x') decays, which we shall 
refer to as the near- field coherence length Axn- In gen- 
eral the result of a measurement of G in this setup is a 
convolution of the object transmission function with the 
near field correlation function Ft^, so that that Axn sets 
the spatial resolution for the reconstruction of the image. 
Note that in the ghost image experiment a bucket detec- 
tion scheme is employed in arm 1, so what we observe in 



,(4) 
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Fig|2t is / dxiG{xi,X2), which makes the imaging inco- 
herent 9] . 

In the ghost diffraction experiment we simply remove 
the lens from the setup of Fig^ so that the detection 
plane of beam 2 is in the focal plane of the lens F. Hence, 
^2(^2,^2) = {iXF)~^e~^^'^'^'^ . From Eq.(^, we obtain 



der^(e,x2)T* 



(-1-0^ 



j d^Tf{lx2) 



(7) 



,(8) 



where f{q) = J ^e-'^-^T{x) is the amplitude of the 
diffraction pattern from the object. r/(x, x') is the sec- 
ond order correlation function ((2|) of the speckle field in 
the far-field plane, as observed in the focal plane of the 
lens F; its correlation length Axf (the far-field coher- 
ence length) sets the limit for the spatial resolution of 
the diffraction pattern reconstruction. 

Relevant to the resolution of the ghost imaging and 
ghost diffraction schemes are hence the spatial coherence 
properties of the speckle field in the near field immedi- 
ately before the object, and in the far-field plane. These 
can be investigated by measuring the fourth-order cor- 
relation functions, in the absence of the object. In the 
near-field plane, we measured the spatial auto-correlation 
function of the reference beam {I2{x)l2{x')) in the setup 
with the lens F' inserted, so that the reference beam 
recorded by the CCD is the (demagnified) image of the 
near-field. This is plotted in Fig^ (squares) as a function 
of \x — x'\. Neglecting the small shot noise contribution 

far-field correlation 
Gaussian fit o=7.8 |im 
■ near-field correlation 
Gaussian fit o=14.3 |im 
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FIG. 4: Fourth order auto-correlation function in the near 
field plane (squares) and in the far- field plane (triangles). The 
full lines are Gaussian fits of the correlation peaks, and the 
data have been normalized to the baseline values. 

at X = x', by using Siegert formula for Gaussian field 
statistics and the BS tranformation 

{h{x)l2{x')) - {I2{X)){I2{X')) 

(4(f)62(fO)f = \r\^ |r,(mf,mfOI^9) 

The baseline in Fig 0] corresponds to the product of mean 
intensities, while the narrow peak 3i x = x' is the r.h.s 



of (13, and reflects the spatial coherence properties of 
the field. By making a Gaussian fit of this peak we 
obtained a variance dn = 14.3 /im. From this we can 
infer the coherence length in the near-field plane, as 
Axn ~ 2mcFn = 34.3 /im. The triangles in Fig^ plot 
the intensity correlation function in the far-field plane, 
obtained by measuring the auto-correlation function of 
beam 1 in the focal plane of the lens F. The nar- 
row peak in this plot is now \t\^ \ V f{x^x')\^ and reflects 
the spatial coherence properties of the speckles in the 
far-field plane. Its variance calculated by means of a 
Gaussian fit is <j/ = 7.8 /im, from which we infer a far- 
field coherence length Ax^ = 2<j/ = 15.6 /im. This in 
turn corresponds to a spread in transverse wave vectors 

I^Ax/ = 1.93 X 10-^ /im-^ 
In the spirit of the EPR argument brought forward by 
[15 Q , we find a product of the near-field and far-field 
resolutions 



AxnAq = 0.066 



(10) 



which is by far smaller than the limit of unity imposed by 
EPR-like arguments in[l3. H^. and is at least four times 
smaller than the results there reported, where entangled 
photons were used to retrieve the ghost image and the 
ghost diffraction pattern. 

Notice that there is no violation of any EPR bound for 
such a separable mixture. In fact, in any plane, the con- 
ditional probability of detecting a photon at position X2 
in beam 2 given that a photon was detected at position 
Xi in beam 1 is: 



^1(5*1)) 



Hence, given that photon 1 was detected at xi there is 
a large probability of detecting a photon anywhere in 
beam 2 (the first term at r.h.s. of Eq. ([TT|) ). with a nar- 
row probability peak localized at X2 = xi. Thanks to the 
first term, there is no non-locality ingredient in the cor- 
relation of the classical speckle beams, as obviously must 
be. However, the crucial point is that the resolution of 
imaging schemes based on correlation measurements de- 
pends only on the widths of the peaks in Fig. ^ and 
that their product in the near and in the far field are not 
bounded by any EPR-like inequality. As a matter of fact 
the near and the far field coherence lengths, correspond- 
ing roughly to the size of the speckles in the near and 
in the far-field, are independent quantities. In the near 
field, the size of the speckle depends on the laser diame- 
ter Dq^ and from the total optical path length z from the 
source to the near-field plane, Axn ex Xz/Dq The 
pinhole D being rather close to the near-field plane, we 
checked that the speckle size in this plane did not depend 
on the pinhole aperture. The pinhole aperture instead in- 
fluences the size of the speckles in the far field, this being 
roughly given by Axf oc XF/D^15\. Using the values 
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of our setup, we find Axn ^ 30 /am and Axf ^ 17 /im, 
in good agreement witii tiie values estimated from tiie 
correlation (Fig. 

In conclusion, we have reported on an experimental 
demostration of high resolution ghost imaging and ghost 
diffraction by using a single source of classically corre- 
lated thermal light. The distributed imaging character of 
our experiment is evident from the fact that the informa- 
tion on the object spatial distribution can be processed 
by only acting on the reference beam 2 (e.g. passing 
from the image to the diffraction pattern). The product 
of resolutions of the ghost image and the ghost diffrac- 
tion schemes is well below the limit that was claimed by 
former literature to be achievable only with entangled 
beams. This could be further improved by optimizing 
the scheme (e.g. with a larger pinhole, a larger laser di- 
ameter) which was constrained in our experiment by the 
need of using a single CCD and a low-power laser. This 
definitely proves the claim set forth in 7], that the only 
advantage of entanglement with respect to classical cor- 
relation lies in the better visibility of information. This 
might be useful in very high sensitivity measurement or 
in quantum information schemes, where e.g. the infor- 
mation need to be hidden to a third party, but does not 
give any evident practical advantage when attempting to 
process information from a macroscopic classical object, 
as in the experiment shown here. 
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(STVF) and the Carlsberg Foundation. 



[1] T. B. Pittman, Y. H. Shih, D. V. Strekalov, and A. V. 

Sergienko, Phys. Rev. A 52, R3429 (1995). 
[2] D. V. Strekalov, A. V. Sergienko, D. N. Klyshko, and 

Y. H. Shih, Phys. Rev. Lett. 74, 3600 (1995). 
[3] P. H. Souto Ribeiro, S. Padua, J. C. Machado da Silva, 

and G. A. Barbosa, Phys. Rev. A 49, 4176 (1994). 
[4] A. F. Abouraddy, B. E. A. Saleh, A. V. Sergienko, and 

M. C. Teich, Phys. Rev. Lett. 87, 123602 (2001). 
[5] A. F. Abouraddy, B. E. A. Saleh, A. V. Sergienko, and 

M. C. Teich, J. Opt. Soc. Amer. B 19, 1174 (2002). 
[6] A. Gatti, E. Brambilla, and L. A. Lugiato, Phys. Rev. 

Lett. 90, 133603 (2003). 
[7] A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, 
Phys. Rev. A 70, 013802 (2004); Phys. Rev. Lett., in 
press; quant-ph/0307187 
[8] M. Bache, E. Brambilla, A. Gatti, and L. A. Lugiato, 

Phys. Rev. A, in pres; quant-ph/0402160 (2004). 
[9] M. Bache, E. Brambilla, A. Gatti, and L. A. Lugiato, In 
preparation (2004). 
[10] E. Brambilla, A. Gatti, M. Bache, and L. A. Lugiato, 

Phys. Rev. A 69, 023802 (2004). 
[11] R. S. Bennink, S. J. Bentley, and R. W. Boyd, Phys. Rev. 

Lett. 89, 113601 (2002). 
[12] J. C. Howell, R. S. Bennink, S. J. Bentley, R. W. Boyd, 

Phys. Rev. Lett. 92, 210403 (2004). 
[13] R. S. Bennink, S. J. Bentley, R. W. Boyd, and J. C. 

Howell, Phys. Rev. Lett. 92, 033601 (2004). 
[14] M. D'Angelo, Y-H Kim, S.P Kulik, and Y. Shih, Phys. 

Rev. Lett. 92, 233601 (2004). 
[15] J. W. Goodman, in Laser speckle and related phenomena, 
edited by J. C. Dainty (Springer, Berlin, 1975), vol. 9 of 
Topics in Applied Physics, p. 9. 
[16] This is only approximately true because the two-lens sys- 
tem is equivalent to a thick lens rather than a thin lens. 



